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Abstract 

As a function of energy E, the average strength function S(E) of a doorway state is 
commonly assumed to be Lorentzian in shape and characterized by two parameters, the 
peak energy Eq and the spreading width Pk The simple picture is modified when the 
density of background states that couple to the doorway state changes significantly in an 
energy interval of size rk For that case we derive an approximate analytical expression 
for S(E). We test our result successfully against numerical simulations. Our result may 
have important implications for shell-model calculations. 
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1. Motivation 

Giant Resonances are an ubiquitous phenomenon in nuclei 0, 0] • A specific nuclear 
mode with normalized wave function |0) carrying definite quantum numbers (spin, parity, 
isospin) is excited, for instance, by absorption of a gamma quantum with specific multi- 
polarity, by nucleon-nucleus scattering, or by stripping of a nucleon from the projectile in 
the collision of two nuclei. The mode has a typical mean excitation energy Eq of several 
or even 10 to 20 MeV, i.e., may occur above the first particle threshold. The Giant Dipole 
(GD) mode in nuclei is a paradigmatic case. Aside from a normalization factor, the wave 
function |0) of the GD mode is the product of the dipole operator and the eigenfunction 
of the nuclear ground state, and the dependence of Eq on mass number A is empirically 
given by Eq « 80 A~ 1 / 3 . In general, the wave function |0) is not an eigenstate of the 
nuclear Hamiltonian H and the mode is, therefore, not observed as a sharp and isolated 
resonance. Rather, the mode spreads in a very short time (typically ps Ti/5 MeV 
w2x 10~ 22 sec) over the eigenstates \i) of H carrying the same quantum numbers (each 
state \i) corresponding to an eigenvalue £i of H). Thus, for the particular reaction under 
consideration the mode |0) acts as a "doorway" to the eigenstates of H which manifests 
itself as a local enhancement of the dependence on energy E of the strength function 

S(£) = £Ko|;>| 2 «5(£- £i ) . (i) 

i 

The levels Si are actually particle-unstable and, thus, resonances, and in most cases S(E) 
is, therefore, a smooth function of E. Often S(E) displays a broad maximum. Pending 
the modifications introduced below, the peak energy is then identified with the mean 
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excitation energy Eg of the doorway state, and the width is identified with and 
referred to as "spreading width" . At excitation energies of f=s 10 MeV, the mean spacing 
d of the nuclear levels e, is typically of order 10 eV, so that ]> 3> d. Hence the name 
giant resonance. Similar phenomena also occur in condensed-matter physics where the 
strength function S(E) is commonly referred to as the local density of states. 

In the simplest theoretical model [3J for the giant-resonance phenomenon, the door- 
way state |0) is coupled to a set of background states \p) (where p = 1,...,N and 
N — > oo) via real coupling matrix elements V^. The background states have constant 
level spacing d. The matrix elements are Gaussian-distributed random variables with 
zero mean values and a common variance v 2 . The strength function is calculated in the 
limit N — > oo as the average over the distribution of the and given by Q 

rV(2*r) 



S(E) ~ (E~E ) 2 + (i/4)(rir ' (2) 

The bar denotes the ensemble average. The average strength function has Lorentzian 
shape and is normalized to unity. The spreading width is given by 

r 4 - = 2nv 2 /d . (3) 

Although it looks like Fermi's golden rule, the result © is correct beyond perturbation 
theory, i.e., for all values of v 2 /d 2 . 

The level density p(E) may be taken to be constant when the rate of change with en- 
ergy of p(E) over an energy interval of length T 4 - is negligible, i.e., when [d\n p(E) / dE]^ 1 <C 
T^. In nuclei, that is not always the case. By way of example we consider the GD mode 
in 16 0. In the shell model |0) is a superposition of one-particle one-hole states. Through 
the residual interaction |0) is coupled to two-particle two-hole states (two particles in 
the sd-shell and two holes in the p— shell). The maximum spacing in energy of the single- 
particle states (of the single-hole states) is about 5 MeV [J] (3 MeV, respectively) , giving 
the spectrum of the two-particle two-hole states a spectral range of about 15 MeV. The 
residual interaction widens the range to ps 25 MeV. The shape of the spectrum being 
Gaussian, the width a of the Gaussian is then around 15 or 20 MeV, and the ratio r 4 «5 
MeV to a is around 1/3 or 1/4 and, thus, not negligible. In the present paper we show 
how Eq. @ is modified under such circumstances. 

Our investigation was triggered by a result for the strength function of a doorway 
state obtained in Ref. [f| . There we considered a Hamiltonian matrix of the form 

The doorway state |0) at energy Eq is coupled to N background states p with p — 
1, . . . , N and N — > oo via real matrix elements V M . The background states are described 
by a real-symmetric random Hamiltonian matrix H^, a member of the Gaussian Or- 
thogonal Ensemble (GOE) of random matrices. The average level density of has the 
shape of a semicircle. Using the Pastur equation we calculated analytically the average 
strength function (the ensemble average of S(E) in Eq. ([T])). Whenever the value of the 
spreading width given by Eq. ([3]) was not negligible in comparison to the radius 2A of 
the GOE semicircle, the effective spreading width r e fi (defined as the full width at half 
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maximum of the average strength function) turned out to be bigger than T^, the increase 
being proportional to T^/X. The method of derivation in Ref. [5[ was confined to the 
GOE with its unrealistic semicircular spectral shape. In the present paper we present an 
approach that, although more approximate than that of Ref. [Hj], applies for a coupling 
of the doorway state to background states with a general dependence of the average level 
density p{E) on energy E. We determine how the effective spreading width T e g differs 
from as given by Eq. ([3]) when p(E) is not constant. 

The model of Ref. [3J disregards all details of nuclear structure. In a more realistic 
approach, one has to replace the statistical assumptions on the matrix elements and 
the assumption of a constant level spacing d by a nuclear-structure model like the shell 
model and/or one of the collective models. In these approaches, the damping mechanism 
has received considerable attention [1, [H 0] , with special focus on the GD resonance 0] • 
Because of the large number of states that couple to the doorway state, the effort is 
substantial, however, and the simple statistical model of Ref. i.e., the use of Eq. @ 
together with a calculation of from Eq. ([3]) , continues to play an important role in the 
analysis of giant -resonance phenomena in nuclei. For that reason we revisit and extend 
the model in the present paper. 



2. Model 

Similarly to Eq. Q we model the doorway state by the Hamiltonian matrix 

E V v 



where the index jj, ranges from 1 to N with N ^> 1. The matrix ^ differs formally from 
that of Eq. ((4]) in that T-L^ has been diagonalized. Instead of the statistical assumptions 
on the matrix H.^ v made below Eq. ((4]), we assume that the are Gaussian random 
variables with zero mean value and a second moment v 2 , and that they are not correlated 
with the E^. We do not need any assumptions on the distribution of the latter. Thus, 
our model is more general than the random-matrix model of Ref. [B| . 
To calculate the strength function S(E), we rewrite Eq. ((T|) as 

S(E) = -l^(0\^-^\Q)) (6) 

where E + = E + ie with e positive infinitesimal. Using Eq. ([5]) we obtain 0, Q 

S ^ = ^(l^^ V(J?+ F WT/ ) ' (7) 

7T \e+ -e -22^^+ -e m ) %y 

Prior to calculating the ensemble average of S(E) we calculate the ensemble average of 
the sum over /i in the denominator of Eq. ((JJ). That sum is denoted by S. The average 
over the distribution of the gives 



E V =v 2 Y I (E+-E,)-K (8) 
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For the remaining sum over p we write 



Averaging over the distribution of the E^, we replace 8{E' — E p ) by p(-E), the average 
level density of the background states, and obtain 

= -tnp(E) + J dE' -JL_p(E') (10) 
where V indicates the principal-value integral. Thus, 

E = -mv 2 p(E) + v 2 J dE' -^L_p{E') . (11) 

We show presently that for N — > oo and 3> d the average strength function S(-E) is 
obtained by replacing in Eq. Q the function E(-E') by E. That yields 

1 T- 1 

S ^ = 2^ (E - E - A) 2 + (l/4)(ri) 2 

where 

= 2nv 2 p(E) , 

A = v 2 J dE' p(E') . (13) 

Eqs. (U2J) and (fT3")) obviously generalize Eqs. © and ([3]) to the case where p(E) is not 
constant and reduce to the latter if it is. 

To justify our averaging procedure (replacement of E by E) we consider first the 
average of S(E) over the Gaussian-distributed matrix elements V p . We use the property 
that the average of the product V^VvVpVa has the value (v 2 ) 2 [5^S pcr + 5 pp 5 vcr + 8 prT S V p], 
and similarly for higher-order terms. In other words, averages over products of Gaussian 
random variables are calculated by Wick contraction of all pairs. For each pair the 
average is zero unless the indices are equal. The average of S(E) in Eq. ([7]) can be 
calculated by expanding the denominator in powers of the Vs and using Wick contraction. 
After averaging, the leading contribution to each term of the series is the one where Vs 
appearing pairwise under the same summation over p are averaged. All other Wick 
contractions restrict the independent summations over p and lead to terms that are 
small of order 1/N and, thus, negligible for N ^> 1. Hence to leading order in 1/N 
averaging S(E) in Eq. Q over the V p is equivalent to averaging E. 

We turn to the average over the E^ and use that S(E) depends on the E p only via 
the expression 5(E — E p ). Expanding S(E) in powers of E and averaging over the E p 
we see that our averaging procedure is justified if the E p are uncorrected. Then, in each 
term of the series J2 P S{E—E p ) is replaced by p(E) and the result is the same as replacing 
E in S(E) by E. The strongest known correlations among eigenvalues are those of the 
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GOE where the follow Wigner-Dyson statistics. GOE level correlations extend over 
an energy range measured in units of d while S(E) varies with energy over an interval 
of length r-k Therefore, such correlations produce correction terms in the expansion of 
S(E) in powers of £ that are small of order d/Y^ and are, thus, negligible for >• d. 
The argument does not apply near the end points of the spectrum where the level density 
tends to zero and d becomes large. This suggests that for our approximation to be valid 
the distance of E from the end points of the spectrum should be larger than T^. We 
observe, however, that equations (|12[) and (|13[) provide reasonable approximations to our 
numerical results even when that condition fails. 

We conclude that Eqs. (|12"j) and (|T3|) for the average strength function S(E) of a 
doorway state are valid except perhaps near the end points of the spectrum. These 
equations generalize Eqs. and ^ by the appearance of a shift function A(E). As 
shown by the second of Eqs. (fT5|) that function accounts for level repulsion between 
the doorway state and the background states. The function A(E) receives negative 
(positive) contributions from background states that lie above (below) the energy E. If 
the spectrum is symmetric about E — then A(0) = and A(E) < (A(E) > 0) if 
E < {E > 0, respectively). For a doorway state at E — this fact widens the spectrum 
and causes r e g to be larger than as given by the first of Eqs. ([15]). Obviously, our 
result agrees with Eqs. @ and © if the average level density of the background states is 
constant so that A = 0. The shift function A(E) is very similar to the shift function for a 
scattering resonance due to its interaction with a continuum of scattering states [^lilfioj]. 

We display the dimensionless ratio A(E)/T^(0) for two important examples: The 
average level density has the shape of a semicircle (the case of the GOE) or of a Gaus- 
sian (this is typical of level densities in the shell-model Q). With the normalization 
/ dE p(E) = N we have 



E ,2 



p(E) = (semicircle), 

p(E) = — ^exp[-£7(2A 2 )] (Gaussian). (14) 
v 2n\ 

Here A denotes half the radius of the semicircle (the variance of the Gaussian, respec- 
tively). With x = E/X the ratio A(E)/T i (0) is given by 

+ 2 



A(£)/r+(o) = -!- / dx' - V— JY^Ta 

Lit J _ 2 x — x' 



I - [9(2 +x)- 6(2 - ~ (semicircle), 



1 r+°° v 

A(£)/r+(0) - — / cLr' ^^exp[-x' 2 /2] 
2tt J_ qo x-x' 



= \e- x2/2 (-i)etf(i^=) (Gaussian). (15) 

Fig. [2 displays these two functions versus x. W conclude that A(E) is significant when- 
ever r-^/A is not negligibly small. 



5 




Figure 1: Plot of the functions A(E)/F^(0) in Eqs. (T5J versus E/X for the semicircle (dashed, color 
online: red) and the Gaussian (solid, color online: black). 



To estimate the effect of the values of A(E) displayed in Fig. [T] on the strength 
function, we note that for the semicircle, the function A(E) is linear in E over the entire 
range of the spectrum while in the Gaussian case, A(E) is approximately linear near the 
center of the spectrum. For both cases we write 

A(E) « axT l {0) = aE (T l {0)/\) = aEj (16) 

where 7 = T^(0)/X and where a > is dimensionless. Substituting that expression into 
Eq. ([12"]) we obtain 

S(E) = ^-r^ ^ ; — (17) 

V ' 2^ 1-07 (E - E ) 2 + V ' 

where 

E - Eq/(1 - orf) and T cS = T l /(1 - aj) . (18) 

The factor 1/(1 — 07) > 1 shifts the mean energy Eq of the doorway state towards smaller 
(larger) values when Eq < (Eq > 0, respectively) and increases the effective value of 
the spreading width. The value of a is obtained by differentiating A(E) at E = 0. For 
the semicircle we find a = 1/4, in agreement with the result of Ref. [Bj]. For the Gaussian 
we have a = l/V2n. In both cases and with 7 w 1/3 or 1/4, that gives a correction of 
about 10 to 20 percent to both Eq and T^. 

In summary, Eqs. (1121) and (|13p are expected to provide a better approximation to 
S(E) than Eqs. @ and ^ if p(E) changes significantly over an energy interval of length 
r-k Then we expect the full width r e g at half maximum of S(E) to be bigger than F-'- 
as given by the first of Eqs. (fl~3|) . Eqs. (|12j) and (|13[) may fail near the end points of the 
spectrum of the background states. This is in accord with the exact results of Ref. Q. 
There it was shown that the interaction with the doorway state increases the range the 
GOE spectrum. Such an effect is beyond the scope of the present approximate treatment. 
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Figure 2: Strength functions (color online: red) of a doorway state at Eq = coupled to background 
states described by a random band matrix of dimension N = 1000 and bandwidth b as defined in the 
text; the coupling matrix elements are uncorrelated Gaussian-distributed random variables (case (i) of 
the text). For comparison, we also display the Lorentzian distribution of Eqs. (£SJ and (color online: 
blue). The average level densities of the background states are also shown (color online: black). Each 
panel corresponds to a fixed value V^(0)/X of the spreading width of Eq. {3}. Averages are performed 
over m = 500 realizations. 



3. Numerical Simulation 

To test the approximations leading to Eqs. (|12p and (|13p we consider a doorway state 
coupled to a random band matrix. Random band matrices have been frequently used 
in different physical contexts. The model is that of Eq. (jU) except that T-L^ v is a real 
symmetric random band matrix of dimension N: All matrix elements with \fi — v\ > b 
vanish. The upper bound on the number of non-zero elements in every row and column 
is (2b — 1). The non- vanishing matrix elements are uncorrelated Gaussian-distributed 
random variables with variances given by Ti 2 ^ = (1 + <5 M! ,)/3 2 . For 6 = 1 the matrix T-L^ v 
is diagonal while for b — N it is equal to the GOE. To make sure that all spectra have 
approximately the same width we determine /3 2 from the condition (l/7V)TraceH 2 = A 2 . 
That gives 

f '^hrY <19) 
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Figure 3: Same as for Fig. [2] but only for 6 = 1 (diagonal random band matrix) and for b = 1000 (GOE 
matrix). In addition to the curves shown previously we now display also the predictions of Eqs. (1126 
and (1131) (color online: blue). 



The average spectrum of the random band matrix is Gaussian for 6=1 and changes 
quickly into an approximately semicircular form as b is increased flU ]. For N — 1000, we 
found that the average spectrum is much more similar to a semicircle than to a Gaussian 
already for b — 5; the transition to semicircular shape was virtually complete at b = 100. 
For b < y/~N and N 3> 1 the eigcnfunctions of a random band matrix are localized, and the 



eigenvalues are uncorrelated, i.e., have Poissonian statistics 11]. Indeed, for N = 1000 
the nearest-neighbor spacing distribution changes from Poisson to Wigner form near 
b = 30. Similarly, the inverse participation ratio defined below decreases strongly with 
increasing b. Some of these results are displayed in the figures shown below. As a 
consequence, random band matrices are useful for testing our approximations both for a 
Gaussian spectrum (b = 1) and for a spectrum with Poisson statistics (b < yN). 

We have considered two ways of coupling the doorway state with the random band 
matrix : (i) The coupling matrix elements in Eq. Q are uncorrelated Gaussian- 
distributed random variables with zero mean values and a common second moment v 2 . 
Then the doorway state is coupled to all states in W^v irrespective of the value of b, i.e., 
irrespective of the localization properties of the eigenvectors of / H. iiV . (ii) We take V M = 
for all [X = 1,,..,N except for /i-values in a band of width w centered in the interval 




Figure 4: Same as for Fig.[2]at Eq/X = — 1, but only for 6=1 (diagonal random band matrix) and for 
b = 1000 (GOE matrix). We also display the predictions of Eqs. il 1 2li and (1131) (color online: blue) and 
of Eqs. Q and J3} (color online: orange). 



[1, N}. The doorway state is coupled only to select states in H^. Localization properties 
of the random band matrix should influence the value of the average strength function 
of the doorway state. In case (ii) the non-vanishing matrix elements were taken to have 
all the same value v^N/w. Then the total coupling strength ^ V% 1S 011 average the 
same in cases (i) and (ii). 

The input parameters of the model are b, w, v 2 , Eq, and N while A defines the spectral 
width and, thus, the energy scale. A further input parameter is m, the number of 
independent drawings of the matrix elements T-L^ v from a random-number generator. 
Each such drawing produces a realization of the random matrix (j4]) . Diagonalization of 
that matrix yields the eigenvalues £j and eigenfunctions \i). These are used to generate 
the strength function in Eq. (TT]). Combining m realizations we obtain the average strength 
function S(E). That function is compared with Eqs. (|12p and ([TU)) . 

For case (ii) with parameter w — 1 the doorway state with energy Eq is coupled to 
only a single other state with Gaussian-distributed energy e, and the average strength 
function S(E) can be calculated analytically. We find 

5(1) = -* ^ _ „, 2 . (20) 



(E - E f 



The function S(E) vanishes at E = Eq, extends over the entire spectrum, and has two 
maxima on opposite sides of E = Eq. That is a consequence of level repulsion and 
explains qualitatively some of the features seen in the figures shown below. 

In Figs. [5] to [5] we present numerical results for case (i). In Fig. [5] we display strength 
functions for parameter values indicated in the figure. The discrepancy between the 
predictions of Eqs. ([2]) and ([3]) and the actual values of the strength function are obvious 
and increase with increasing values of the spreading width T^. We note the gradual 
development of a dip at E = and of a double-hump of the strength function. We 
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Figure 5: Full width at half maximum r e g of the average strength function of the doorway state, 
displayed as a function of the input width IN- (0) of Eq. (|3} ■ The circles represent results of the numerical 
diagonalization, the lines correspond to Eqs. i l 1 21 ) and (1 1 3 D - 

believe that these features correspond to properties of the simple model of Eq. (P20"l) . 
In Fig. [3] we compare some of these results with the predictions of Eqs. (JT2J) and (fT5)) 
and find very good agreement. We note that the dip is correctly reproduced. We have 
performed similar calculations for non-zero values of Eq and found that the Lorentzian 
model shows even larger discrepancies, since it is not able to reproduce not only the 
correct width, but also the asymmetric shape that develops for Eq ^ 0. On the other 
hand, Eqs. (fT2"]l and (flU)) provide the same good agreement for every value of Eq. An 
example is shown in Fig. [4] for the case of Eq/X = —1, while our results are summarized 
in Fig. [5j We believe that the agreement of the numerical results with Eqs. (fT2]l and (fT3|) 
is impressive. We also note that r e ff and differ significantly. 

We turn to case (ii). Again we have performed calculations for N = 1000. Case 
(ii) agrees with case (i) for b = 1000. The results for b = 10 are qualitatively similar 
to those of case (i). Therefore, we focus attention on the case 6=1 and consider 
Eq — 0. Results are shown and compared with the exact analytical result of Eq. ((20)) in 
Fig. [5J The double hump is clearly displayed. The agreement is very good as expected. 
Similarly good agreement was found when Eq was chosen different from zero. On the 
other hand, a comparison of Fig. [5] and the results of Eqs. (fT2"j) and (fTS)) displayed in 
Fig.[3]shows that the approximations leading to Eqs. (fT2"j) and (fT5)l fail when the doorway 
state is coupled to a single background state only (the case w = 1 in Fig [6J . That is 
a very special situation and not typical for doorway states. Increasing the number of 
background states to which the doorway is coupled but keeping the bandwidth of H^v 
unchanged, very quickly changes the strength function so that approximate agreement 



10 




Figure 6: The average strength function (color online: red) for a doorway state coupled to a random 
band matrix with b = 1 of dimension N = 1000 (case (ii) of the text) is compared with the analytical 
result of Eq. (1 20 II (color online: blue). The average level density of the backgound states is also shown 
(color online: black). The value of T^(p)/X is given in each panel. The dashed lines (color online: green) 
show the strength function when the number of states directly coupled to the doorway state is increased 
from one to three. 

with Eqs. (|12D and (fT5| is attained. For w = 3 that is shown by the dashed lines (color 
online: green) in Fig. [6] 

To understand the role of localization in the mixing of the doorway state with the 
background states we have calculated the average inverse participation ratio (IPR) of the 
doorway state. The IPR is defined in terms of the amplitudes (0\i) of the expansion of 
the doorway state in the basis of eigenfunctions \i), i = 0, 1, . . . , N of the matrix (0]) as 
J2i I (0|i) I 4 - If the doorway state is spread more or less uniformly over the eigenstates then 
the normalization condition |(0|i)| 2 = 1 suggests that the IPR has a value around 
1/N. If, on the other hand, the doorway state mixes with only a few of the eigenstates 
then the IPR should be much larger than 1/N. Thus, for case (i) we expect values 
of the IPR around 1/N and for case (ii) much bigger values. The left panel of Fig. [7] 
corresponds to case (i). The IPR (black dots) decreases as increases. The values for 
6 = 1 are somewhat larger than those for b = 1000 but still close to 1/N. The solid, 
dashed and dashed-dotted lines are obtained from the simple estimate d/T^(0) for the 
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Figure 7: Inverse participation ratio of the doorway state as a function of the input spreading width 
r^/A for N = 1000 and several values of b as shown in the panels. In the left panel the coupling matrix 
elements are uncorrelated Gaussian-distributed random variables (case (i) of the text), while in the right 
panel they are equal and restricted to a band of width w = b (case (ii) of the text). The average level 
spacing of each matrix ensemble is denoted by d. 

IPR. For case (ii) (right panel) and 6=1 and 6 = 10 the IPR is significantly larger than 
1/N and roughly given by 1/6. That shows that the doorway state mixes only with few 
(ps 6) states. 

4. Conclusions 

We have investigated the strength function of a doorway state coupled to a number 
of background states in cases of strong coupling (the spreading width is not small 
compared to the range of the spectrum of background states). Our result in Eqs. (TT21 and 
(|13p generalizes the standard weak-coupling result and agrees with it for weak coupling. 
We have tested our result by numerical simulations. In most cases studied, we found 
perfect agreement between the numerical results and Eqs. (fT2"| and (fT5|) . Exceptions are 
found only when the doorway is coupled to a single background state. That situation 
is atypical. Even when the number of directly coupled states is increased from one to 
three, approximate agreement with Eqs. (1121) and (I13[) is attained. 

We have pointed out that the strong-coupling case is of practical interest and may 
actually play a role in shell-model calculations. That is the case whenever the spreading 
width is not very small compared to the range in energy over which the average level 
density of the background states changes significantly. Then our Eqs. p2[) and (fT3"|) offer 
a more accurate description of the strength function of the doorway state than do the 
standard Eqs. © and (J3J). Typically the full width at half maximum r c g- of the average 
strength function is then larger than the theoretical expression = 2irv 2 p(E). The 
difference may be important for a comparison between theory and data. 

When the doorway state couples to a random band matrix with localized eigenfunc- 
tions and when the coupling involves only a narrow band of states (our case (ii)), we have 
found that the inverse participation ratio of the doorway state is considerably larger than 
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the inverse matrix dimension. That shows that the doorway state mixes only with a re- 
stricted number of localized states. The result is important for practical applications. 
Indeed, in the nuclear shell model the doorway state (a lp lh state) mixes with 2p 2h 
states which have a Gaussian spectrum. But because of the presence of other modes of 
excitation, the actual nuclear spectrum is not Gaussian in shape but increases monotoni- 
cally with energy, and one may ask what significance our results have in view of this fact. 
However, the mixing of the 2p 2h states with such other states is weak (otherwise shell 
structure would not persist). Modeling such weak mixing in terms of a random band 
matrix with localization, we have shown that our results remain valid in the presence of 
other modes of excitation. 
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